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' When the underlying stock price is a strict local martingale process under an 

^ ■ equivalent local martingale measure, Black-Scholes PDE associated with an European 

option may have multiple solutions. In this paper, we study an approximation for 
O^l the smallest hedging price of such an European option. Our results show that a 

class of rebate barrier options can be used for this approximation. Among of them, 
, a specific rebate option is also provided with a continuous rebate function, which 

I corresponds to the unique classical solution of the associated parabolic PDE. Such a 

00 ' construction makes existing numerical PDE techniques applicable for its computation. 

, An asymptotic convergence rate is also studied when the knocked-out barrier moves 

' to infinity under suitable conditions. 

CM 
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1 Introduction 

In a financial market equipped with the unique equivalent local martingale measure (ELMM) 
P, the smallest hedging price of an European option is the conditional expectation of the 
payoff with respect to the probability P, see [HI [7] . In contrast to the probabilistic repre- 
sentation, option price can be also characterized as the unique solution of its associated 
Black-Scholes PDE, provided that PDE has a unique classical solution. 

The necessary and sufficient condition for the unique solvability of the parabolic PDE is 
that, the underlying stock price is Martingale process, see [1]. In other words, if the stock 
price is a strict local Martingale, then there exists multiple solutions for Black-Scholes PDE. 
Moreover, the option price may be one of the many solutions, see |1] . The difference of the 
multiple solutions of PDE is termed as financial bubbles, see [21 HI [6] and the references 
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therein. In this work, we will consider the following problem proposed by Fernholz and 
Karatzas 0: 

(Q) How can one find a feasible numerical solution convergent to the option price under 
the lack of uniqueness of Black- Scholes PDF? 

We first examine the existing numerical schemes on CEV model of Example [U where the 
option price can be explicitly identified. There are typically two kinds of numerical schemes 
in this vein [11] . One is Monte Carlo method by discretizing the probability representation, 
the other is PDE numerical method by discretizing the truncated version of PDE. 

Unfortunately, Example |2] and Example [3] shows that classical Euler-Maruyama ap- 
proximation (for Monte Carlo method) and finite difference method (for PDE numerical 
method) leads to a strictly larger value than the desired option price. Motivated from these 
two examples, question (Q) boils down into the following two problems: 

(Ql) Find a feasible approximation for a Monte Carlo method, and its convergence rate; 

(Q2) Find a feasible approximation for PDE numerical method, and its convergence rate. 

In short, this work intends to find a feasible approximation to the smallest superhedging 
price V{x,t). It turns out that the value function can be obtained by a limit of a series of 
appropriate rebate option prices, which can be estimated by usual Monte Carlo method, see 
the details in Corollary El However, Corollary [2] may not be utilized for the approximation 
by PDE numerical method, since it may cause a discontinuity at the corner of the terminial- 
boundary datum. Therfore, a specific rebate option is proposed with its price continuous 
up to the boundary, so that its price corresponds to the unique classical solution of its 
associated parabolic PDE. Such a construction makes existing numerical PDE techniques 
applicable for computations, see Theorem |3l 

The rest of the paper is outlined as follows. In the next section, we give precise formu- 
lation of the problem. Section 3 presents main results, and related proofs is relegated to 
Section 4 for the reader's convenience. The last section summarizes the work. 

2 Problem formulation 

Throughout this paper, we use K as a generic constant, and M"^ = (0, oo), M"*" = M+ U {0}. 
If A is a subset of M x [0, T], then C{A) denotes the set of all continuous real functions on 
A, C'^'^{A) denotes a collection of all functions ip : A M. such that ip^x and ipt belong to 
C{A). D^{A) denotes the set of all measurable functions ip : A ^ M"*" satisfying growth 
condition 

ip{x,t) < K{l + \x\^), 'i{x,t) e A. (2.1) 

C^{A) = C{A) r\D^{A) denotes the set of all continuous functions satisfying 7-growth. We 
also denote the parabolic domain Q := M"*" x (0, T), truncated domain := (0, /3) x (0, T), 
and := (a, /3) x (0, T) for < a < /3. 
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We consider a single stock in the presence of the unique equivalent local martingale 
measure (ELMM) P, under which the deflated price process follows 

dX{s) = a{X{s))dW{s), X{t) = x > 0, (2.2) 

where is a standard Brownian motion with respect to a given probability space (^2, J-", P, F : 
{^s '■ s >t}) satisfying usual conditions. We impose the following two conditions on / and 
a: 

(Al) a is locally Holder continuous with exponent | satisfying a{x) > for all x G M.~^, 
a{0) = 0. 

(A2) / : 1+ ^ 1+ is a C^(t+) payoff function for some 7 G [0, 1]. 

By [HI 5.5.11], the assumption (Al) on a ensures there exists a unique strong solution of 
( 12. 2 p with absorbing state at zero. 

For a contingent claim f{X{T)) with a fixed maturity T > 0, the smallest hedging price 
has the form of 

Vix,t) = E,.,,[/(X(T))] := E[f{X^'\T))\Tt]- (2.3) 

In the above, we suppress the superscripts {x,t) in X^'*, and write Kx.t[ ■ ] to indicate the 
expectation with respect to P computed under these initial conditions. 

Recently, [1] shows that the value function V of 02.31) is the C'^'^{Q) fl C{Q) solution of 
BS{Q,f), where BS{Q,f) refers to Black-Scholes equation 

f (E) Ut + ^a\x)u,, = on g = M+ X (0, T) 

[ {BD) M(0,t) = /(O) on Mt G (0,T]. 

However, the next example taken from [2] shows that the value function V may not 
be the unique solution of BS{Q, f) when the deflated price process X is a strict local 
martingale. 

Example 1 (CEV model). Suppose the stock price follows a strict local martingale process 
dX{s) = X^dW{s), with the initial X{t) = x > 0. Consider V{x,t) = E^.t[X(r)]. Then, 
V can be computed explicitly as 

V(.,t)=.{l-2<,{--±=)). (2.5) 

One can verify V satisfies BS{Q, /). Another trivial solution is u{x,t) = x. 

Now, V is one of the possibly multiple solutions of BS{Q,f). With the existence of 
multiple solutions to PDE BS{Q, /), our question is 

• If we use any of existing PDE numerical methods on BS{Q,f), or any of existing 
Monte-carlo methods on (12. 3p . does it converge to the desired value function among 
multiple solutions of PDE? 
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Unfortunately, the answer is NO in general. In fact, the next trivial example shows that 
the classical Monte Carlo method by Euler-Maruyama approximation does not lead to the 
desired value V{x,t) of (12 .Sp of Example [H 

Example 2. Consider the strong Euler-Maruyama (EM) approximation to Example^with 
step size A, 

^n+i = + a{X^)iWinA + A) - W{nA)), = x. 

Let X^(-) be the piecewise constant interpolation of {X^ : n > 0}, i.e. 

^^(^) = ^[t/A], Vs>0. (2.6) 

Since {X^ : n > 0} is a martingale, the approximated value function simply leads to a 
wrong value 

VA{x,t) := E,,t[X^{T)] = E,,4X(^_,)/J =x> V{x,t). □ 

Similar to Monte-carlo method. One can also prove that the finite difference method on 
PDE BS{Q, /) also leads to a wrong value. 

Example 3. Black-Scholes PDE associated to the CEV model Example^ is BS{Q,f) of 
(12. 4p with a{x) = x^ and f{x) = x. To use the finite difference method (EDM) in the 
above PDE, we shall truncate the domain and put artificial boundary conditions on the 
upper barrier {{x,t) : <t < T} for large enough /3 > 0. ^45 suggested by /Ti] /. we impose 
boundary conditions, which asymptotes the option price, i.e. 

u{^,t) = /3, VO < t < T. 

With step size in space variable x and A in time variable t, one can easily check upward 
finite difference scheme backward in time yields trivial numerical solution u^{x,t) = x for 
any small A > 0. □ 

To the end, our work is to resolve the following question: How can one find a feasible 
approximation of this value function V of (12.31) in both Monte Carlo method and EDM? 
What is the convergence rate? 

3 Main result 

In this subsection, we present the main results, and the proofs will be relegated to the next 
section. 
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3.1 Approximation by Monte Carlo method 

We consider the following up-rebate option prices: Suppose the up barrier is given by a 
positive constant f3 > x > and stopping time t'^ (suppressing the initial condition {x,t)) 
is the first hitting time of the stock price X{s) to the barrier /3, i.e. 

^x,t,/3 _ ^ ^ . j^x,t^^^ >/3}AT. (3.1) 

For some function g, let its payoff at r^'*''' consist of 

1. rebate payoff g{f3), if < T; 

2. otherwise, terminal payoff f{X{T)). 
Then, the rebate option price is of the form 

V^{x,t) = E,M(3n{r^<T} + f{X{T))l{^,^T}], (3.2) 

is a functional of / and g, and we may write V^'^'^ instead of whenever it needs 
an explicit emphasis on its dependence of / and g. It turns out that the the option price 
value of (12. 3 p can be obtained by a limit of a series of appropriate rebate option prices. 

Theorem 1. Assume (A1-A2). Suppose the rebate payoff g satisfies one of the following 
two conditions: 

1. g{x) is of sub-linear growth, i.e. lim^j^oo = 0; 

2. g{x) is of linear growth, i.e. limsup^_^go < ^> ^'^'^ ^■^ ^ martingale. 
Then, we have the convergence for of (13. 2p . 

lim V^{x,t) = V{x,t). 

In addition, if g G /^^(M"*") with'-fArj < 1, then the convergence rate is the order o/l — (7X777) 
as (3 ^ 00, i.e. 

I (V - V^) {x, t) I < A'/J-^^-^^^"" , Vx < /3. (3.3) 

Theorem [1] shows that Monte Carlo method on the expression Vf^ of (Q actually 
leads to correct estimation of the option price V, provided that the rebate function g is 
appropriately chosen. Among the many choices of g, the simplest one shall be taking (7 = 0. 
Next result summarizes the above comments. 

Corollary 2. Let 

y'^'° = E.,,[/(X(T))l|,,=r]. 

Then, limp^ooV^'^lXjt) = V{x,t) point wisely in x and t. Furthermore, if f{x) = 0{x'^ 
with some constant 7 < 1, then its convergence rate is 

\{V - V(^){x,t)\ = 0{p-^+^), as (3 ^00. 
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In this below, we fix Monte Carlo method in Example [2] based on convergence result in 
Theorem [1] 

Example 4. Let's extend f : R to f : R ^ R by f{x) = /(O) for x < 0. With 

^'^i') of (12.61) . let be the first hitting time of to the barrier (5. The modified Monte 
Carlo scheme to approximate V{x^t) of Example\^is given by, 

V^{x,t) :=E,,i 

Corollary^ implies that the rebate option price is convergent to the smallest hedging price, 
i.e. 

V^{x, t) := E,,i[/(X(T))l|,,>r}] ^ V{x, t) as (3 ^ oo. 

Note that 

X^iT) ^ X{T) a.s. and \f{X^{T))l , \ < max /(x). 

I A— J (J<X<p 

Hence, Bounded Convergence Theorem also implies that 

V^ix, t) y/^ix, t), as A ^ 0. (3.4) 

As a result, the option price V{x,t) of (12. 3p can be approximated by using Monte Carlo 
method on the above rebate options, i.e. 

lim lim V^f(x,t) = V{x,t). 

/3— >oo A->-0 

Regarding the estimation by Monte Carlo method, one may take the simplest choice 
g(f3) = for the rebate payoff as of Corollary [2l see also [3]. However, Corollary [2] can not 
be utilized the approximation by PDE numerical method, since it may cause a discontinuity 
at the corner (/3,T) of the terminial-boundary datum when /(/3) 7^ 0. 



/(X^(T))1^^. 



3.2 Approximation by PDE numerical method 

For the above Monte Carlo method on (13.21) . yet another to be mentioned is a drawback 
in the computation by PDE numerical methods due to the possible discontinuity of the 
boundary-terminal data. 

To illustrate this issue, we write Black-Scholes PDE associated to the rebate option 
price V^ix^t) of flO) . 

Ut + ^a\x)u,, = on := (0,/3) x (0,T); 

M(x,T) = /(x), VxG[0,/3]; (3-5) 
M(0,t) = /(0), u{(3,t)=g{^), VtG(0,r). 

Note that, PDE (13. 5p has a discontinuous corner at the point (/3,T) if g{(3) 7^ f{/3)- Also 
recall that, the choice oi g{P) = /(/3) may not be possible, like in CEV model of Example [1] 
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It is well known that, if g{f3) ^ /(/?) and the boundary-terminal data is discontinu- 
ous, then one can not expect the unique solution of (13.51) continuous up to the boundary. 
Furthermore, the discontinuity and the singularity at the corner propagate the numerical 
errors quickly throughout its entire domain for the numerical PDE methods, such as finite 
element method (FEM) or finite difference method (FDM), see more discussions in [10] and 
the references therein. Therefore, the unique solvability and the regularity of the solution 
are crucial to make use of the existing PDE numerical methods. 

To avoid this error propagation due to the discontinuity of the boundary-terminal data, 
we provide an alternative choice to (13. 2p by revising the terminal payoff: Consider a rebate 
option of barrier /3 with 

1. zero rebate payoff, i.e. g{f3) = 0; 

2. and a revised terminal payoff 

fix) = /(x)l{.<,/2} + ^^^^^f^lw2<.m- (3.6) 

In this case, the rebate option price 

V^{x,t) = E,,4/(X(T))l{,,=r|] (3.7) 

is associated to PDE 



(^)/3 Ut + ^a\x)u,, = on := (0, /3) x (0, T); 

(BD) m(0, t) = /(O), uiP, t) = 0, Vt e (0, T); 
{TD)^ u{x,T) = f{x), Va;e [0,/3]. 



(3.^ 



Observe that, the revised terminal data /'^ not only makes the terminal-boundary data 
continuous at the corner (/3, T), but also preserves Holder regularity of the original terminal 
data / regardless how large the value /3 is. 

Although PDE (13. 8 p is degenerate at x = 0, one can still has unique classical solution 
by utilizing Shauder's interior estimate. Also, its solution is indeed equal to the revised 
rebate option price of (13. 7p . see Lemma HI Moreover, by using comparison principle 
twice on two different truncated domains, one can show its unique solution V'^ must be 
convergent to the desired value V of (12. 3p . 

Theorem 3. Assume (A1-A2). Then, of (jST]) is the unique C'^'^{Qi3)nC(Q/^) solution 
of PDE ([31]), and 

lim V^{x,t) = V{x,t), V(x,t) e Q. 
In addition, if j < 1 in (A2), then the convergence rate is 

-V\{x,t) < Kp-^+\ 

Thanks to the Theorem [3l one can use either well established FDM or FEM on PDE 
(13. 8p for a large (3 to estimate the smallest superhedging price. 
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4 Proof of main results 



In this section, we will first characterize the value function V. Based on the properties 
of we can estimate \V — V^l to prove Theorem [1], and \V — V'^l to prove Theorem [3l 
respectively. 

4.1 Characterization of the option price V 

We have seen that the option price V of (12.31) is one of the solutions of BS{Q,f). To 
proceed, we need identify which solution corresponds to the option price V among many. 
This enables us to establish the connection between parabolic partial differential equation 
BS{Q, f) and probability representation (12. 3p . 

Proposition 4. Assume (A1-A2). Then, value function V of (12.31) is 

1. the smallest lower-hounded C'^'^{Q) fl C^{Q) solution of BS{Q, /). 

2. the unique C'^'^{Q) fl C{Q) solution of BS{Q, f) if and only if a satisfies 



Proof Theorem 3.2 of [1] shows that Visa C^^\Q)nC{Q) solution of BS{Q, /). Applying 
super-martingale property of X(T) and Jensen's inequality, the next derivation shows that 
V G C^(Q), 



For the necessary and sufficient condition on uniqueness, we refer the proof to pj. It 
remains to show V is the smallest lower bounded solution. Sometimes, we use X to denote 
X^'" without ambiguity in this proof. Note that, by path- wise uniqueness of the solution 



to (O 

Y{t) ^ v{x^'\t),t) = E[/(x^^'°W'*(r))|j-d = E[/(x-'°(r))|j-,] 

is a martingale process. Suppose V G C'^'^{Q)nC{Q) is an arbitrary lower bounded solution 
of BS{Q, /), then Ito's formula applying to Y{t) = V{X{t),t) leads to 




(4.1) 



V{x,t) = E,,i[/(X(T))] < A^l +E,,i[X^(T)]) < K(l + xT). 



y(t) = r(X(0),0)+ / V,{X{s),s)aiX{s))dWis), 



and Y{t) is a lower bounded local martingale, hence is a super-martingale. Therefore, we 
have 



Y{0) > E[F(T)] = E[/(X(T))] = r(0) 



and this implies 



V{x,0) > V{x,0). 

We can similarly prove for V{x,t) > V{x,t) for all t. □ 
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Proposition 5. Assuming (A1-A2), BS{Q,f) only admits non-negative solution in the 
space of lower bounded C'^'^{Q) fl C{Q) functions. 

Proof. Proposition H] shows that V is the smallest lower-bounded solution of PDE. Since 
V > hj definition of f l2.3p . it implies any lower-bounded solution u satisfies u > V > 
0. □ 

In Example [H we have seen that BS{Q, f ) of CEV model has multiple solutions. We 
continue this model to demonstrate Proposition 0, a solution smaller than V must be 
unbounded from below. 

Example 5. By Proposition^ the explicit solution V >Q of (12. 5 p in CEV model smallest 
lower-hounded solution of BS{Q, /). In fact one can find, 

is a smaller solution, i.e. v <V in Q. However, v is not lower-bounded, i.e. v{x,t) — ?■ — oo 
as X — 7- oo. 

4.2 Proof of Theorem [1] 

Recall that the domain of the value function V is given on the domain Q of (12. 4p . and its 
related truncated domain is given by (13. 5p . Let (f : Q ^ be a measurable function. 
We introduce the truncated value function V^'"^ for convenience, 

[ f{x,t) Otherwise. 

where the stopping time of (13. ip is the first hitting time to the barrier (3. By the above 
definition, 

= V^{x,t) 

for the of (E^D, if we set 

ipi{x,t) = g{x)l[t<T} + /(x)l{f=T}- (4.3) 

With the above setup, to prove Theorem^ our goal is to estimate \ V^''^^ — V\ as (3 ^ oo 
with ifi of (14. 3p and the constraint on g given in Theorem [H We emphasize here, ipi may 
not be continuous up to the boundary, i.e. ipi ^ C{Q) when g{x) < f{x) for some x > 0. 

Lemma 6. Assume (A1-A2). Then, 

1. V{x, t) = V^'^ix, t) for alio < X < (3. 

2. If if , ijj : Q —> M"*" are two measurable functions satisfying ip > tp on d*Qp, then 
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Proof. X^'* is the unique strong solution of (12. 2p due to (Al). Therefore, the conclu- 
sion follows from the following simple derivation using tower property and strong Markov 
property: 

V{x,t) =E[fiX^'\T))\T,] 

= E[/(X-'*(T))|J-.,]|J-,] 

= E[E[/(X^(-')'-'(r))|J-,,]|J-,] 

= E[ViXiT^),T^)\T,] 

= V^^^{x,t). 

Monotonicity of V^'"^ in Lf follows directly from the definition of V^'"-^ of (14. 2p . □ 

It is noted that, two results of Lemma [6] correspond to uniqueness and comparison 
principle of its associated PDE. However, we provide the probabilistic proof Lemma [6] , 
since we want to cover potentially discontinuous function ipi of (14. 3p . in which uniqueness 
may not remain true. 

Lemma 7. Assume (A1-A2) andg>0. Then, V^^^^ defined by (g^]) and (jOD satisfies 

lim V^^'^'ix^t) > V{x,t). 



In addition, equality holds in the above if and only if 

lim E^,t g{l3)l^rt<<T} 

/3^-oo L 







(4.4) 



Proof. We start with the following observation: The solution X := X*'^ of (12. 2 p does not 
explode almost surely by [HI 5.5.3], i.e. 

lim = T, a.s.-P (4.5) 

/3->oo 

Due to this fact together with Monotone Convergence Theorem, we obtain following iden- 
tities: 



lim/3^oo ^x,t 



f{X{T))^^,^r} 



\imp^ocf{X{T))1^^0^T} 
'f{X{T))]=V{x,t). 



(4.6) 



By the definition of ipi of (14. 3p . this results in 
lim V"^''^i(x,t) 

= lim E,,t[^,iXiT^),T^)l{,,^T}]+ lim E,,i[v9i(X(rO, r^l^^^^^}] 

/3^cx> /3-i-oo 

= lim E,,i[(^i(/3,r^)l|,,<r}] + lim E,.,i[/(X(T))1{,,=^}] 
= lim E^4gi(3)l^,p^T}] + Vix,t). 

Rearranging the above identity, we have 

V{x,t) = lim V^''^^{x,t) - lim E,,,[^7(/3)1{.^<t}]- 

/3— >-CXD ^— 5>oo 



(4.7) 



Note that three terms in (14. 7p are all non-negative. Hence, limj3^ooV^''^^{x,t) > V{x,t) 
and equality holds if and only if (14. 4p holds. □ 
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As mentioned in fl4.5p . the solution X^'* of (12. 2p does not explode almost surely, and 
this can be rewritten as 

p(^^,t,/3 < T) ^ as /3 ^ oo. 

An interesting question about this is that, how fast does the above probability converge to 
zero? The answer to this question is indeed useful to obtain the convergence rate of the 
truncated approximation. 

Proposition 8. Fix (x, t) G Q and assume (A1-A2). As fi ^ oc, stopping time r^^'^'i^ of 
(13. ip satisfies 

1. P{r^'*'/3 <T} = 0(l//3). 

2. Moreover, P{r^'*'^ < T} = o{l/(3) if and only if {X^'^{s) : t < s <T} is a martingale 

Proof. By taking g{x) = f{x) = x in (14.71) . 

lim E,,i[X(r^)] = lim /3P{r=^'*'/^ < T} + E,. i[X(T)]. 

For all P > X, since {X^'*(r^ A s) : s > t} is a bounded local martingale, hence it is 
martingale. So, E2,.^f [X(r^)] = x for all /3 > x. Rearranging the above identity, we have 

lim /3P{r^'*'/' < T} = x - E^,i[X(T)] (4.8) 



( I48|) implies 

1. Since E^t[X{T)] > 0, lim^j^oo < T} < x < oo, which shows P{r^''*'^ < T} = 
0(1/(3): 

2. {X*'^'(s) : t < s < T} is a martingale if and only if x = E^ t[X(T)], if and only if 

P{r^'*'^ <T} = o(l//3). 

□ 

Finally, we are now ready fo the proof of Theorem [H 
Proof of Theorem [1]. We first show its convergence, then obtain convergence rate. 
1. Regarding its convergence, it is enough to verify (14. 4p by Lemma [3 Note that 



lim E^.^t 



gm^^.^T}] < lim ^ lim (3F{t^'''^ < T}. 
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(a) If g is of sub-linear growth, then lim/3^00 = 0. Hence, (14. 4p holds due to the 
first result of Proposition El 
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(b) On the other hand, if X*'^ is a martingale, then we have hm^^oo /3P{r^'*''' < 
T} = from the second result of Proposition [HI and fl4.4p remains true provided 
that g is of linear growth. 

2. Since V{x^t) = V^'^{x,t) for all /3 > x by Lemma [6], we have the following identity: 

Setting V := sup^g^.T) ^i^^t): we can rewrite 

\{V-V^nM\ < {\V\ + \g\){mMr^<T}]. (4.9) 

Since V G C^(M+) by Proposition H and g E D^{R+), we have \V + g\ G D^y^{R+). 
Hence, write f l4.9p by Proposition |8] 

\{V - V^'^^){x,t)\ < (|F| + \g\m\0{l/^) < 

which finally results in f l3.3p . 

□ 

4.3 Proof of Theorem [3] 

Lemma 9. Assume (A1-A2). Then of f l4.2p is the unique solution of (13. 8p in the space 
ofC''\Q^)nC(Q^). 

Proof. Fix G Qi^. Take a G (0,x/2). Recall = <5/3 H {Q^Y be an open set. Also 

define 

r"'^ = inf{s>t:(X^'*(s),s)^g^}. 
Due to the uniform ellipticity, 

V"'^ix,t):=E,,[f^iX{T'^'^))] (4.10) 

is the unique classical solution of 

ut + la\x)u^^ = 0, on Q° = (a, /3) x (0, T) 

w(/3, t) = 0, M(a, t) = f^ia), Vt G (0, T) (4.11) 
«(a;,T) = /^(x), Va;G[a,/3]. 

If we restrict V^'^ on the subdomain Q^J'^, it solves following PDE uniquely, 

Ut + |(t2(x)m^^ = 0, on g^/^ = (x/2, /3) x (0, T) 

u{P, t) = 0, m(x/2, t) = V'''^{x/2, t), Vt G (0, T) (4.12) 
m(x,T) = /^(x), Vx G [x/2,/3]. 
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Furthermore, by using Shauder estimate Theorem 4.9 together with Theorem 5.9 of [9], 
one can have estimate on weighted Holder norm, i.e. 

for some constant K independent to a. On the other hand, by definition f l4.10p . we have 

/2 = supV^"''^ < sup \f^{x)\ < sup \f'^{x)\<K 

for some K independent to a. Let d = \ minja;, /9 — x, t, T — t}, which must be less than 
the minimum distance of x to any point in the parabolic boundary d*Q'^. Consider a 
neighborhood of x given by Nx{r) := {x — r,x + r) x (t — r, t + r). By definition of the 
weighted norm (Page 47 of [9]), we finally have the following a-uniform estimate on Nx{d), 

Therefore, Arzela-Ascoli Theorem implies that there exists a subsequence of : a G 

{0,x/2)}, which is uniformly convergent to a function u on Nx{d), i.e. 

V^'^ — J- u as a — )■ 0, uniformly on Nx{d). 

The uniform convergence implies that the limit function is u G C'^'^{Nr^{d)). Using the facts 
of almost sure convergence r"'^ — ?• r^, together with dominated convergence theorem, one 
can check that 

lim V'^ix^t) = E^,41im fiXir""'^))] = V^{x,t), pointwisely. 

a— s>0 ' Q— >-0 

Hence, = u solves {E)p of (13. 8p in the classical sense. By bounded convergence theorem, 
one can also show V^{x,t) G C{Qp) from the facts 

lim VHx, t) = 0, lim V^(x, t) = /(O), lim V^(x, t) = fix), (4.13) 

x^P a;->0 t^T 

Thus, we conclude is the classical solution of (13. 8p . Moreover, strong solution satisfies 
maximum principle, and hence the uniqueness follows from Corollary 2.4 of [5]. □ 

Now, we are ready to prove Theorem [31 

Proof of Theorem [3l V^^ is the unique solution of (13. 8 p by Lemma [9l Fix (xo,to) ^ Q 
and /3 > 2xo. We will use comparison principle of Lemma [6] twice to obtain the desired 
results. Define (y92 : Q t-> M by 

^2{x,t) = f{x)l{t=T}- 
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Since ip2{x,t) < on d*Qi3/2, we can apply Lemma[6]on to obtain V^^'^''^^{xo,to) < 
V^{x(),tQ). Similarly, since < (f2{x,t) on d*Qi3 by its definition, we apply Lemma on 
Qi^ to obtain V^''^^{xo,tQ) > V^{xo,to). Thus, we have inequality 

V^/'''''{xo,to) < V^{xo,to) < V^''^'{xo,to). (4.14) 

Taking lim/3_=.oo in the above inequality and using Theorem [1], all three terms shall converge 
to the same value V{xo,tQ). The rate of the convergence is the combined result of (14.141) 
and dajD- □ 



5 Further remarks 

This paper studies an approximation to the smallest hedging price of European option using 
rebate options. From mathematical point of view, this work concerns on the approximation 
of the value function V of (12. 3p by truncating the domain Q and imposing suitable Cauchy- 
Dirichlet data g. 

The main result on the convergence Theorem [T] provides that, if the function g is chosen 
to satisfy sublinear growth in x uniformly in t G [0,T), then the truncated value V^'^ 
converges to V. This enables practitioners to adopt EM methods on big enough truncated 
domain Q/? to get a close value of V^, as demonstrated in Example HI 

On the other hand, to adopt numerical PDE techniques, continuous Cauchy-Dirichlet 
data is desired to get a good approximation. However, if the payoff / is given as of a 
linear growth, g is taken as of a sublinear growth in x for the purpose of the convergence 
by Theorem [H then it's not possible to have a continuous solution of Black-Scholes PDE. 
Alternatively, we provide a continuous Cauchy-Dirichlet data by modifying the terminal 
payoff appropriately. 
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